Nonequilibrium Kondo Effect in a Quantum Dot Coupled to Ferromagnetic Leads 
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We study the Kondo effect in the electron transport through a quantum dot coupled to ferro- 
magnetic leads, using a real-time diagrammatic technique which provides a systematic description 
of the nonequilibrium dynamics of a system with strong local electron correlations. We evaluate the 
theory in an extension of the 'resonant tunneling approximation', introduced earlier, by introducing 
the self-energy of the off-diagonal component of the reduced propagator in spin space. In this way 
we develop a charge and spin conserving approximation that accounts not only for Kondo correla- 
tions but also for the spin splitting and spin accumulation out of equilibrium. We show that the 
Kondo resonances, split by the applied bias voltage, may be spin polarized. A left-right asymmetry 
in the coupling strength and/or spin polarization of the electrodes significantly affects both the 
spin accumulation and the weight of the split Kondo resonances out of equilibrium. The effects are 
observable in the nonlinear differential conductance. We also discuss the influence of decoherence 
on the Kondo resonance in the frame of the real-time formulation. 

PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 



I. INTRODUCTION 

The continuing experimental progress with mesoscopic 
electronic devices has stimulated anew the interest in fun- 
damental quantum mechanical questions A good 
example is the Kondo effect in electron transport 
through quantum dots (QD). The effect, which nowa- 
days is well established experimentally 0, has remained 
a challenging problem for nonequilibrium quantum trans- 
port theory (ill EH- 

After early experiments, which 
concentrated on semiconductor QDs, the Kondo effect 
was observed in single-atom Q and single-molecule tran- 
sistors , as well as in carbon nanotube QDs ^3 j all of 
which were coupled to metallic leads. More recently, the 
Kondo effect was demonstrated in QDs coupled to ferro- 
magnetic leads ^3 ^|. In this article, we will discuss 
the nonequilibrium Kondo effect in such systems as an 
example of the nontrivial collective many-body physics 
in magnetic nanostructures. 

Research on electronic properties in magnetic systems 
has evolved into an active field, "spintronics" 13]. In 
spintronics devices magnetic properties control transport 
properties via the electron spin degree of freedom. One 
example is the tunnel magnetoresistance (TMR) in fer- 
romagnetic tunnel junctions [l3 |. which relies on switch- 
ing the magnetization of the electrodes between paral- 
lel (P) and antiparallel (AP) orientations by an applied 
magnetic field. In magnetic mul tilay ers, the intrinsic ex- 
change coupling between layers [ISj influences the mag- 
netic structure and consequently the transport proper- 
ties . The Kondo effect in transport through a QD 
coupled to ferromagnetic leads (Fig. ^ provides a prime 
example of many-body effects in spintronics devices. The 
theoretical model, a local spin coupled to two ferromag- 
netic leads, was originally introduced to explain the de- 



crease of the TMR with increasing temperature ■ The 
model was revived recently, and led to a series of publica- 
tions m El m El E il mi . in the early stage, due 
to different approximations used, the question remained 
controversial whether the Kondo resonance splits in the 
presence of spin-polarized leads [H E3 '^o* IH 113 ■ 
The splitting was confirmed eventually by numerical- 
renormalization-group (NRG) techniques [2a. l2^ . 

The origin of the splitting of the Kondo resonances lies 
an exchange coupling induced by spin-dependent quan- 
tum charge fluctuations. Let us consider a single- level, 
originally spin-degenerate QD with energy eg < (mea- 
sured relative to the Fermi level of the lead electrons) 
and strong Coulomb interaction, U oo. Then, dou- 
ble occupancy is suppressed and the QD behaves as a 
magnetic impurity. But there still exist quantum charge 
fluctuations due to the tunneling between the QD and 
the leads. An electron with majority-spin in the QD 
can tunnel between the QD and the leads easier than an 
electron with minority-spin, thus gaining a larger kinetic 
energy. As a result, a ferromagnetic exchange interac- 
tion between the QD spin and the lead magnetization 
is induced. It lifts the spin degeneracy of the QD level 
by Ae. Such a spin splitting quenches spin-flip scatter- 
ing processes for low energies, below |Ae|, and suppresses 
the Kondo correlations. It also leads to a splitting of the 
Kondo resonance. The qualitative two-stage scaling the- 
ory presented in Ref. [13] indicates that the high-energy 
part of quantum charge fluctuations, between |eo| and a 
cutoff energy D (of the order of the itinerant electron 
band width) is responsible for a spin splitting given by 
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Here we assumed the coupling strength Po- between the 
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QD level and the leads to be spin dependent. The split- 
ting of the level results in a corresponding spin splitting 
of the Kondo peak in the spectral density. Remarkably, 
it is possible to compensate the spin splitting Ae by ap- 
plying an external magnetic field, and thus to recover the 
full Kondo effect 

In spite of the interest in the problem, a systematic 
approximation which can account for the spin splitting 
of the Kondo effect in nonequilibrium transport problems 
has not yet been presented. 

- The NRG is a powerful tool for the investigation of 
the Kondo problem, but is limited to the equilibrium 
state. The validity of extensions to the nonequilibrium 
state ,!£, 2Sl 22^ remains to be proven. 

- The slave-boson mean-field approximation is correct in 
the limit of large degeneracy and able to capture the 
physics of the Kondo singlet at zero temperature [25l |. 
However, it does not account for the spin splitting since 
the description of quantum charge fluctuations is beyond 
mean- field theory. 

- The non-crossing approximation (NCA), a systematic 
conserving approximation ^^G^ ,27j , has been successfully 
adopted to the nonequilibrium state 0,01 . NCA accounts 
for spin-dependent quantum charge fluctuations and thus 
for the spin splitting due to ferromagnetic leads. But, it 
also generates an additional spurious peak at the Fermi 
level even in equilibrium pol | as well as a spurious peak in 
the presence of magnetic field . Such a shortcoming is 
dangerous in the context of transport, which is sensitive 
to the electron states near the Fermi energy. The spuri- 
ous peak can be made to disappear, if vertex corrections 
are properly accounted for |28j| . However, the general- 
ization of NCA described in Ref. does not appear 
tractable in nonequilibrium situations. 

- The equation-of-motion (EOM) approach, which was 
adopted in Ref. 0, is able to capture Kondo correla- 
tions as well as - in principle - the spin splitting. How- 
ever, within the standard decoupling scheme, the spin 
splitting is lost 18, 21], and it requires introducing an ad- 
ditional self-consistent equation to determine the renor- 
malized QD-level energy [l^ . 

In this situation, a systematic approximation, free of 
the mentioned drawbacks in nonequilibrium, is needed 
for further progress in the exploration of the Kondo ef- 
fect in spintronics devices. In this paper, we formulate 
the problem using a real-time diagrammatic technique, 
which is a systematic method suitable to describe the 
nonequilibrium time evolution of a system with strong lo- 
cal electron correlations f^, '2^ . We evaluate the theory, 
accounting for Kondo correlations and spin-dependent 
quantum charge fluctuations, in a conserving approxima- 
tion making use of an extension of the resonant tunneling 
approximation (RTA) J]. The RTA was shown to pro- 
duce qualitatively reasonable results in the presence of 
a magnetic field at not too low temperatures, sufficient 
to allow for a comparison with experiment [3y|. This 
property motivated us to apply the RTA for the present 
problem. 



(a) 



t 



(b) 



R 



t t 



R 



FIG. 1: A single level quantum dot coupled to two ferromag- 
netic leads for (a) parallel and (b) antiparallel alignments of 
lead magnetizations. We allow for a left-right asymmetry in 
tunnel barriers and/or spin polarization factors. 



The outline of this paper is as follows. In Sec. IHl we 
introduce the model Hamiltonian and shortly review the 
real-time diagrammatic technique as well as the resonant 
tunneling approximation. We introduce an extension of 
the RTA, which is needed to account for the spin split- 
ting. In Sec. nni we present numerical results. We will 

show the splitting of the zero-bias anomaly for various 
temperatures fSec. IhT^ . We provide a comprehensive 
discussion of the restoration of the Kondo resonance by 
an applied magnetic field f Sec, mi's)) and the nonequilib- 
rium Kondo effect fSec. IIII H)! . In Sec. IIII Dl we will dis- 
cuss the effect of asymmetries in system parameters and 
their consequences for nonequilibrium states. In Sec. lIVI 
we will discuss decoherence effects on the Kondo reso- 
nance out of equilibrium within the framework of the 
real-time diagrammatic technique. Sec. [V] is devoted to 
a discussion of the relation between our results and recent 
experiments ^3^3- We summarize in Sec. IVII 



II. FORMULATION 
A. Model Hamiltonian 

Our model Hamiltonian consists of several parts de- 
scribing a single-level QD, Hu, the left (and right) reser- 
voir. 



and tunneling, Ht, 



H = Ht 



E 

■■=L.,R 



Hr + 



(2) 



The Hamiltonians for the isolated QD and the leads are 
Hu=Y^ eadU^ + Ud\d^d\di , (3) 

(T 

Hr — ^ Srka alka^^ka , (4) 
rk(T 

respectively (we use unit where h — kB — I)- Here, 
da is the annihilation operator of electrons with spin-cr 
(cr^t, I) in the QD, and U denotes the onsite Coulomb 
interaction. The energy of the single QD level with spin 
t/i is ey i = Co ± Ez/2 + aceV/2. The first term, eo, 
may be tuned by a gate voltage. The Zeeman energy is 
given by Ez = 2giiBB, where the external magnetic field, 
B, is assumed to act (only) on the spin in the QD. The 
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applied bias voltage between the left and right leads, V, 
may shift the QD level by aceV/2 Here is an asym- 
metry factor due to different capacitive couplings of the 
QD to the left and right leads. The operator Urka annihi- 
lates an electron in the lead r — L,R with wave number k 
and spin a. We assumed that the ferromagnetism in the 
leads can be treated within mean-field approximation, 
i.e., it is accounted for by a spin-dependent dispersion 
£rka- and, hence, spin-dependent density of states (DOS) 

The tunneling Hamiltonian, finally, is given by 

Ht^YI ^rk dlarka + H.c. (5) 

rk(T 

We assumed that the tunneling amplitude Vrk is indepen- 
dent of spin. Furthermore, for simplicity, we will ignore 
a k dependence of the tunneling amplitude Vrk = Vr. 

The important parameter of this model is the spin- 
dependent coupling strength to the ferromagnetic lead r, 

r,,(Lj) = 27r|K-|'!^r.(w) . (6) 

The spin dependence is the main generalization required 
for the present problem as compared to that of a QD 
coupled to normal metal leads 7]. The parameter r^, 
which appeared in Eq. is given by To- = J^r^ra, 
where Tra = ^ra-iO) is the value at the Fermi energy. The 
spin-dependent coupling strengths are related to the spin 
polarization of the leads by 

-(r,t-r,j/(r,| + r,i). (7) 

In this paper, we will analyze two configurations of the 
lead magnetizations, the parallel (P) and antiparallel 
(AP) alignments. We always take a positive value for 
Pl- The P (AP) alignment then corresponds to positive 
(negative) values of Pr. 

B. Real-time diagrammatic technique 

For infinitely strong Coulomb interaction, U = 00, dou- 
ble occupancy of the QD level is suppressed, and the 
Hilbert space is restricted to the empty state |0) and two 
singly occupied states \a). In the presence of strong in- 
teraction effect, the application of Wick's theorem is pro- 
hibited, and the standard approach for the calculation of 
noncquilibrium currents through a tunnel junction [sH 
is not suitable. We, therefore, proceed using the 'real- 
time diagrammatic technique' developed by Schoeller et 
al. 29]. Apart from providing a transparent classifica- 
tion of various processes it allows constructing a conserv- 
ing approximation out of equilibrium. In the following, 
we will briefly outline the technique. It is closely re- 
lated the Keldysh formalism and the Feynmann- Vernon 
approach (see for example Ref. [33). The technique en- 
ables one to perform a systematic diagrammatic expan- 
sion of the reduced density matrix in terms of iJx in the 
real-time domain. 



The central quantity is the reduced density matrix, 
which is derived from the total density matrix by tracing 
out the lead electron degrees of freedom, 

P^, (t) = Tr^Tr^Trz, [p^p^p? Tc^c^x' : («) 

X^,^(i)7^e'^^«*|x')(x|e-*^"*, (x,x' = 0,T,i). (9) 

The combination of a trace over dot states, Tr^, and a 
projection operator, Xy^> ^, selects the (x, x')"Component 
of the reduced density matrix P^,. The subscript / de- 
notes the interaction picture with respect to the unper- 
turbed Hamiltonian Hq = Hl + Hj^ + Hu- The contour 
ordering operator Tc is defined on the closed time-path 
C starting from the initial time, ti, to time t on the up- 
per branch of the Keldysh contour and returning to ti on 
the lower branch. The Keldysh S-matrix is defined along 
this contour by 

Sc = exp (^-i J dtHT{t)i^ . (10) 

Initially, at time U, the density matrix is assumed to be 
factorized into a product of equilibrium density matrices 
of the leads, pg, and the density matrix of the QD, pp . 
The former are characterized by the temperature T and 
chemical potential p^ ifJ-L = —fJ-R = £V/2) as, 

p5 = e-''(^'--S'=-^-°'^-"'-^"V^J , (11) 

where P = T~^. The normalization factor Zq ensures that 
Tr,.p[; = l. 

Equation (jHJ can be evaluated in a perturbative ex- 
pansion in terms of Ht- Under the conditions that the 
initial density matrix of the QD is diagonal, (x'|Po*lx) = 
^x'^x'.X' ^'^^ that the Hamiltonian conserves the spin, 
the time evolution reduces in the limit of — cx3 to a 
master equation: 

E ^x'^x.x-O, (12) 

x=o,hl 

for the (stationary) probability P^ which is the diago- 
nal (x, x)"Component of the stationary density matrix. 
The self-energy of the master equation S^.x' related to 
transition probabilities from a state |x) to a state |x')- 
The main task is to determine an appropriate self-energy 
of the master equation in the frame of the diagrammatic 
expansion. 

Figure |21 shows examples of diagrams. Basic build- 
ing blocks are the free evolution (horizontal lines on the 
Keldysh contour) and the tunneling processes (remaining 
directed lines) of electrons with spin a into or out of the 
QD to the leads r 0. In Fourier space the latter carry 
a factor obtained by Fermi's golden rule, 

7± (^) = rrA^)f^{uj - p.)/(2^), (13) 

where f^{uj) = l/(e*^'^-|-l) is the electron (hole) distribu- 
tion function. In oder to regularize the expressions we im- 
pose a Lorentzian cutoff replacing the coupling strengths 
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FIG. 2: Diagrams for various tunneling processes: (i) quan- 
tum charge fluctuations, (ii) the co-tunneling process and (iii) 
the resonant tunneling process. The horizontal lines on the 
Keldysh contour represent the free propagators of the quan- 
tum dot. Remaining directed lines are tunneling lines. 



of Eq. © by TraPc{uj) where pduj) = / {uj^ + D^) . The 
effective bandwidth D was already introduced in Eq. Q . 

In Fig. 121 each sector of the diagram corresponds to 
a physical tunneling process. The 'bubbles' in sector (i) 
represent quantum charge fluctuations, where an electron 
with spin a virtually fluctuates between the QD and the 
leads. A single line connecting the upper and lower hori- 
zontal lines (not shown) describes a sequential tunneling 
process. Two overlapping parallel lines shown in sector 
(ii) represent a co-tunneling process, where in one coher- 
ent transition an electron with spin a leaves the QD and 
an electron with spin a' fills the empty level. A 'train' of 
overlapping lines [sector (iii)] represents a process with 
electrons going back and forth many times between the 
QD and the leads. 



C. Extension of resonant tunneling approximation 

In order to account for the logarithmic behavior of 
Kondo correlations, we have to sum up a suitable class of 
diagrams of the expansion up to infinite order. The start- 
ing point of our approximation is the class of diagrams 
denoted as 'resonant tunneling approximation' (RTA) 
In the diagrammatic language, RTA includes all diagrams 
where any auxiliary vertical line cuts at most two tun- 
neling lines [An example is presented in Fig. [21 iii)]. At 
not too low temperatures RTA was shown to provide 
a qualitatively reasonable approximation for the Kondo 
correlations with analytical expressions for the nonlin- 
ear differential conductance for arbitrary level degener- 
acy. It predicts also qualitatively reasonable behavior in 
the presence of Zeeman splitting. However, the origi- 
nal RTA developed in Ref. does not describe the spin 
splitting since, as we will discuss below, it does not take 
into account the subclass of diagrams which accounts for 
spin-dependent quantum charge fluctuations. 

To find the origin of the shortcoming, let us examine 
the QD spin dynamics. In general, in order to describe 
the superposition of spin states 1 1) and | i), it is necessary 
to consider also off-diagonal components of the density 
matrix in spin space, e.g., , where a is the opposite 
spin of a. The time evolution of between time ti and 
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FIG. 3: (a) Dyson equation for 7r?|^ and (b) the diagram- 
matic expansion of the irreducible self-energy cr^^J. Directed 
dotted lines carry the energy ui. 



t {t > ti) is described by the reduced propagator 

T^s,l{t,ti) = -i{^\TrLTTR [pQPQTcScXg^{t)i] |cr). 

(14) 

Since the tunneling Hamiltonian Eq. |(SJ) conserves the 
spin, spin is conserved at each vertex. E.g., the spin in- 
dices of an incoming/outgoing tunneling line and an out- 
going/incoming free propagator for each vertex are the 
same. Therefore, after Fourier transformation, Eq. H14|) 
can be expressed as a Dyson equation, 

^?;?(c.) = l/ico + iO+ -e, + e^- fT?;^(c.)), (15) 

where (Tg]% denotes an irreducible self-energy and 0^ is 
a positive infinitesimal number. The diagrammatic rep- 
resentation of Eq. H15|l is given in Fig. OJa). Directed 
dotted lines carrying the energy lu correspond to the fac- 
tor exp{iLut) of the Fourier transformation. The diagrams 
for cr^'^ are shown in Fig. Ol^b), in which any auxiliary 
vertical line cuts at least one tunneling line. 

Now, we observe that within RTA in any diagram for 
the self-energy of the master equation Eq. (|12|l a sector, 
where the reduced propagator ttJ^^ appears, is always ac- 
companied by two tunneling lines, a left- going tunneling 
line with spin a and a right-going tunneling line with spin 
a, which are disconnected from tt?'^ itself. It is a con- 
sequence of the requirement that any auxiliary vertical 
line must cut the same number of left- and right-going 
tunneling lines or free propagators with the same spin 
indices. This in turn follows from the fact that the spin 
is conserved at each vertex, and that the initial QD den- 
sity matrix is assumed to be diagonal. As a result, the 
irreducible self-energy does not appear in RTA 0, 
because in its presence, any auxiliary vertical line laying 
on a sector of (J%'% would cut at least three tunneling lines. 
RTA can account for the Zeeman splitting of the Kondo 
resonance, which is already included in a denominator of 
the bare reduced propagator 7r^|^. However, RTA omits 
further information on the spin dynamics. As we will 
discuss in the following, the further information, such as 
the spin splitting Ae caused by spin-dependent quantum 
charge fluctuations, is included in beyond the origi- 
nal RTA. 

In order to obtain an expression for a%'\^ we proceed 
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in a perturbative expansion in terms of the coupling F, 

(16) 



G.a a.(T (1) I (7,a (2) 



In lowest order we find the first two diagrams in the 
right-hand side of Fig. |2Ib) describing spin-dependent 
quantum charge fluctuations. From them we obtain the 
first-order term 



''duo' 



= l-{uJ + eg)-%{e,-^)\ (17) 

where 7^ = 7ro- ^"^^ if — Ylr Irh- Here we introduced 
7^, which is the Hilbert transform of 7^, 



^p(a;-M.){±^Q+2^ 



T Reih{ - + I — 



and '0(x) denotes the digamma function. One can expect 
that the lowest order is sufficient for the description of 
the spin splitting. Indeed, for Ez = T = eV = and 
ferromagnetic leads {Pl—Pb. = P>0), the spin splitting 
Ae obtained by the scaling theory [Eq. J^l] is reproduced, 

crj'l (cj) ~ Ae for Iw/eol <C 1. We further note that the 
self-energy given by Eq. (|17|l also describes the reduction 
of the Kondo resonance splitting below the value of the 
Zeeman energy for normal (nonmagnetic) leads. For 
Pl=Pr = Q and T — eV = the self-energy then reads. 



{lj) k, tanh 

TT 



a.r 



(19) 



for Cct < and Iw/co-l <C 1, which renormalizes and re- 
duces Ez- Such a reduction was derived before for the 
Anderson [ssj and Kondo [s^l models. 

Next, we discuss the nature and details of the extended 
RTA. The extension amounts to inserting into all possible 
positions in the diagrams of the original RTA expansion 
the lowest-order irreducible self-energy a^'^'J'^^ given by 
Eq. (|17|l . The reduced propagator describing the time 
evolution of P^ is given by Eq. H14|) with the replacement 
of a by 0: 



(20) 



The resulting self-energy of the reduced propagator is 
represented by the diagrams from Fig. ^a). It yields 



a. (7 ( \ 



E 

s = ± 



duj'- 



uj + iO+ - uj' 
+ ldiv'^ti.^')^lf'\^-^')- 



(21) 



Here t^%%^^^ is defined by Eq. (|15|l with the self-energy 
given by Eq. (|17|l . The result is similar to the self-energy 
of the Green function obtained previously using EOM 
(Eq. (8) in Ref. (ijI). The first term of Eq. (EU is re- 
lated with spin a [the first and the second diagrams in 
the right-hand side of Fig. El^a)]. The real part of this 
term causes a slight shift of the QD-lcvel and the imag- 
inary part gives a level broadening. The second term 
[the third diagram in the right-hand side of Fig. Ela)] 
is associated with the Fermi distribution function of the 
opposite spin a through 7^ and is roughly proportional 
to / dw/+(w)/(w-FjO+ tEzT Ae) for cr=T/i at eV = 0. 
Consequently, the second term leads to the logarithmic 
behavior of Kondo correlations. Our result for a?'^ pro- 
vides a microscopic justification to lowest-order expan- 
sion in terms of F for the spin splitting crAe postulated 
in the EOM scheme (Eq. (9) of Ref. [13). 

The self-energy of the master equation H12|) is given by 



2zlm f duo J2 A^l^l{Lo)^^^ico) , (22) 

s=± 



and Ex,o = -Jla^i.i ^x,<t (x = 0. T, i)- The vertex func- 
tions A+ and A~ are determined by solving the diagram- 
matic equations depicted in Fig. ^b). In the diagram 
for A+/~, an open line carrying energy uj (dotted line) 
connects to the upper/lower horizontal line 35[. The 
diagrammatic equations can generate diagrams for the 
resonant tunneling processes in a recursive manner. The 
explicit forms are 



s=± 



w + «0+ - uj' 



(T ( /\+ s{ l\ cr.cr(l)/ /\ I 



(23) 



Again in the right-hand side of Eq. (|23|l , we inserted the 
self-energy a^y via TT^y ' . 

From P^', which is obtained by solving the master 
equation Eq. (|12|l , we find the average occupation N and 
the local magnetization M of the QD, 



N^Pf^Pf, M= (l/2)(Pf -Pf) 



(25) 



The spin-dependent nonequilibrium DOS of the QD is 
defined as 



P.H = -^{G>(^)-G<(^)}, 



(26) 



where the lesser/greater Green functions are expressed, 
using Eqs. and 1231), as 

G|(c.) = 2zlm^ i'^'A±^;o(c^). (27) 
x=o,T,i 
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FIG. 4: (a) The irreducible self-energy of the reduced prop- 
agator ttq'q . (b) Diagrammatic equations for the vertex func- 
tions A+^^o a-nd A~iJ^o. An open dotted line connects to the 
upper/lower horizontal line for A^. Shaded boxes represent 
the off-diagonal component of the reduced propagator in the 

o-,<t(1) 

spm space vr^'^ • 



III. RESULTS AND DISCUSSIONS 

In this section, we will present results obtained by 
solving the integral equations numerically. In all results 
we tested the numerical accuracy by checking the spin 
and charge conservation ( the relative error is smaller 
than 5 X 10~^). We also checked the spectral sum rule, 
/ duj pa{uj)=PQ^ + P^*", which follows from the fact that 
A^^\q {lj) is analytic in the upper half plane and decays 
at infinity as 1/uj. The sum rule is always satisfied with 
accuracy better than 5 x lO"**. 

In Sec. nil Al we will discuss the splitting of the zero- 
bias anomaly. Subsections IIII Bl and IIII CI provide a dis- 
cussion of the restoration of the Kondo resonance as well 
as of properties of the nonequilibrium Kondo effect. In 
theses subsections, we will assume a left-right symmetry 
in the coupling strengths Tl=Tr [Fr = (Fr|-f rr^)/2] and 
spin polarization factors Pl = \ Pr \ , and ac — Q. In equilib- 
rium these symmetries are not important because, as one 
can see from a unitary transformation 23] , an asymmetry 
gives just a constant prefactor to the linear conductance. 
However, this is not the case for a nonequilibrium state, 
as will be shown in Sec. IIII Dl 



Also the transport current flowing into the QD through 
the junction r, 7^, can be expressed in terms of the Green 
functions. The spin- resolved components read 

Ir. = -le J doj {7+ (c^) G>{lo) + 7,,- (c.) G<(c.)}. (28) 

One can prove that the approximation satisfies spin 
and charge conservation [s^. Equations (|22|l . H27|l and 
combine to J2r ^ra = « e P^* S^^cr , which by virtue 
of the master equation, Eq. (|12|l leads to the spin-resolved 
current conservation, '^j, Ira = 0. The conservation laws 
for charge and spin, 

Il = -Ir = I, ^(/.t-^h) = 0, (29) 

r 

follows the spin-resolved current conservation. Equation 
(|29|1 is valid only for collinear magnetizations, where the 
spin is conserved at each vertex in the diagrams of the 
self-energy for the master equation, regardless of the sign 
of Pl(r)- This fact ensures that the self-energy of master 
equation does not mix diagonal and off-diagonal states 
and allows one to formulate a closed equation for the 
diagonal components of the reduced density matrix. 

The linear conductance, Gjj" = Yivav^o dIa[V)/dV 
[In =/rcr], probcs the equilibrium spin-dependent DOS, 
P'^[lij), around the Fermi energy: 

where Gk = e^/(27rfi,) is the conductance quantum. We 
symmetrized the expression Eq. (|30|l using the spin and 
charge conservation, Eq. H29|) . 



A. Splitting of the zero bias anomaly 

Figure Efa) shows the nonlinear differential conduc- 
tance Gp{V) =dI{V)/dV for P alignment and symmet- 
ric coupling, Tl — Fij, and spin polarization, P^ = Pr- 
The splitting of the zero-bias anomaly is well observable 
for low temperatures. The enhancement of the peak at 
low temperatures is a clear sign of Kondo correlations. 
The value of the splitting is independent of temperature 
and almost coincides with the prediction of the scaling 
theory Eq. (UJ, 2Aew0.204F. At high temperature, the 
two peaks are smeared out. However, the magnitude of 
the splitting is not affected by the temperature. 

Figure |Sfb) shows the corresponding plot for the local 
magnetization M. For small bias voltage \eV\ < |Ae|, it 
increases with decreasing temperature. At the same time, 
the differential conductance decreases [Fig. (Sfa)]. This 
indicates that the suppression of the zero-bias anomaly 
can be attributed to the quenching of the spin-flip scat- 
tering at low temperatures. When the applied bias volt- 
age exceeds the spin-splitting energy, \eV\ > |Ae|, a spin- 
flip channel opens. As a result, the local magnetization 
starts to decrease and peaks appear in the nonlinear dif- 
ferential conductance. 



B. Restoration of Kondo resonance 

A remarkable feature of the spin splitting induced by 
ferromagnetism is that it can be compensated by an ap- 
plied magnetic fleld To', ^23] . At the point of full compen- 
sation, the strong coupling limit can be retrieved. The 
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FIG. 5: (a) The nonlinear differential conductance and (b) 
the local magnetization M for the P alignment [Pl ~ Pr ~ 
0.1) and for symmetric coupling rL=F[i. Solid, dashed, dot- 
ted and dot-dashed lines are results for various temperatures 
T/r = 0.005,0.01,0.025, and 0.05, respectively. Parameters 
are: eo/r = -2, D/r = 50. 

Kondo temperature in the leading logarithmic approxi- 
mation is given by 

„ , f — 7re arctan(P) 1 
TK{P)^iexp^— J^j- (31) 

Here we considered the case Pj^ — — P. The ratio 
r/(7re) corresponds the value Jo(i^| + i^x) of Ref. 0. The 
cutoff energy e follows from e = — eo + r/(27r) ln(e/£'). 

Figure Efa) and (b) show the linear conductance G'" 
and the local magnetization M as a function of Zeeman 
energy Ei in equihbrium for the parallel alignment. The 
average occupation N depends slightly on Ei and varies 
between iV«0.88 and A^«0.91. At -Ez^O, the direction 
of the local magnetization M and lead magnetizations are 
aligned. With increasing Zeeman energy the local magne- 
tization M changes the direction. At a peak in the linear 
conductance G''", where spin-flip processes are expected 
to be enhanced [solid line in Fig.E^a)], we observe the lo- 
cal niagnetization M = in agreement with results from 
Ref. i23| . The spin state of the QD is responsible for the 
slight asymmetry in the conductance around the M — 
point. In the left (right) side of the peak, because the 
QD level is occupied by an up-spin (down-spin) electron, 
the co-tunneling current for up-spin (down-spin) compo- 
nent is dominant [dashed (dotted) line in Fig. |Hfa)]. As 
the up-spins are majority spins, a larger conductance is 
obtained on the left side of the peak. 

We analyzed the DOS dX E-z=0 [Fig.EIc)] and at Ez ^ 
0.205 F, where Af « [Fig.EJd)]. The solid and dashed 
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FIG. 6: (a) The linear conductance G''" and the local mag- 
netization M as a function of the Zeeman energy Ez for the P 
alignment {Pl = Pr = 0.2) . In panel (a) , the spin-resolved con- 
ductance is plotted by dashed (dotted) line for up (down) spin. 
The equilibrium DOS for (c) Ez/T^O and (d) Ez/r= 0.205. 
Solid and dashed (dotted) lines show the total DOS and the 
spin resolved DOS for up (down) spin. Insets in these pan- 
els are blow-ups of the Kondo resonance. Parameters are: 
T/r = 0.005, eo/r = -2, D/r = 50. 



(dotted) lines show the total DOS p{uj) = J^a P<y{^) 
the up-spin (down-spin) component of the DOS p„(uj), 
respectively. Figure |SJc) and its inset demonstrate that 
the Kondo resonance splits aX Ez — 0- The two res- 
onances, which appear rather well resolved, should be 
smeared by decoherence effects. These effects are be- 
yond the present approximation, as we will discuss in 
Sec. II VI The enhancement /suppression of the up/down 
spin (dashed/dotted line) component of the QD-level 
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DOS (broad peak around cj/r = — 2) indicates the posi- 
tive local magnetization, M>0. 

At Ez = 0.205 r, the Kendo resonance is restored at 
the Fermi energy [solid line in Fig. EJd)]. We can ob- 
serve that the resonance for the minority-spin component 
(dotted line) is higher than that of the the majority-spin 
component (dashed line) [inset in Fig. El^d)]. Loosely 
speaking, it is because the minority (majority) spin in 
the QD is well (poorly) screened by majority (minority) 
spin in leads. This tendency is consistent with the zero- 
temperature result of NRG ,23J and the Friedel sum rule, 
which predicts Po-(O) ~ l/F^ [23|. 

As for the QD-level DOS, up-spin and down-spin com- 
ponents are almost the same (at Ez = 0.205 F), which 
indicates M = [dotted and dashed lines in Figs. EJ^d)]. 
However, though the QD level can be spin polarized de- 
pending on Ez, the total QD-level DOS is not affected 
by Ez, i.e. the solid lines in Figs. EIc) and (d) can al- 
most overlap around the broad peak. Such a feature is 
also obtained from the NRG work [2H. 



C. Kondo resonance out of equilibrium 

In this subsection, we will discuss the nonequilibrium 
Kondo effect. Fig.Qshows the nonlinear differential con- 
ductance, G [(a) and (c)], and the local magnetization, 
M [(b) and (d)], for symmetric coupling. As discussed 
in Sec. IIII Al we can observe the splitting of the zero- 
bias anomaly for parallel alignment [solid line in (a)]. 
However, for antiparallel alignment [solid line in (c)], the 
splitting vanishes. Moreover, for positive (negative) bias 
voltage spin accumulates in the QD with positive (nega- 
tive) local magnetization [solid line in (d)]. 

As discussed in Sec. lHIBl the zero-bias anomaly is ex- 
pected to be restored when the spin splitting is compen- 
sated by an applied magnetic field. The dashed line in 
(c) shows the case for the P alignment with Ez =0.205 F. 
Figures[3(e) and (f) show the nonequilibrium DOS Pcr(w) 
for the AP alignment and that for the P alignment with 
the Zeeman energy Ez = 0.205 F. The total DOS has 
two peaks (solid lines) at two chemical potentials, simi- 
lar in shape to that of a nonequilibrium Kondo DOS split 
by an applied voltage. However, now the two peaks are 
spin polarized, i.e. there is clear difference in the weights 
for the spin-up component Pi{i-o) (dashed lines) and the 
spin-down component pj,(w) (dot-dashed lines). 

For the AP alignment with the Zeeman energy Ez = 
0.205 F, the splitting persists [dashed line in (a)], but in 
this case the splitting is caused by the Zeeman energy ^ , 
and a strong asymmetry between two peaks is obtained. 
Intuitively, the result can be understood by the following 
argument: In the regime \eV\ < Ez, a spin-down occu- 
pies the QD level. Above a positive bias voltage eV>Ez, 
spin-flip tunneling events from a left lead spin-up state 
to a right lead spin-down state contribute to the cur- 
rent. In this case, because both initial and final states 
are majority-spin states, a large peak is expected. On 




- -0.2 



-0.4 -0.2 0.2 0.4 -04 -0.2 0.2 0.4 

ev/r eV/r 



0.2 


1 

AP 






(e) 


p 




1 


1 


(f) 






fy Pf+P-- 














=^0.1 




— • \ KT 

\p. - 















Ez= 0.205r 








Ez= 


o.iosr 








0.2 


H \ 1— 


-1— 




— h 
{g)_ 


1 1 1 

>^ A 


■ 


— h- 


h 

(h)_ 


■=^0.1 


~- ,/\ ,f 














eV=-Ez 







1 1 




1 1 











-0.4 0.4 -0.4 0.4 

co/r co/F 



FIG. 7: The nonlinear differential conductance G, (a) and (c), 
and the local magnetization M, (b) and (d), for a symmetric 
coupling r_t =r_R. Solid lines and dashed lines in (a) and (c) 
are plots for Ez=Q and i5z = 0.205 F, respectively. The solid 
lines in (a) and (c) are for the P alignment [Pl = Pr ~ 0.2) 
and the AP alignment {Pl = ~Pr = 0.2)], respectively. The 
dashed lines in (a) and (c) are for the AP alignment and the P 
alignment, respectively. The average occupation is N ~ 0.89 
for all parameter regimes. The panels (e) and (f) show the 
corresponding DOS p{i^) of the solid and the dashed lines 
in panel (c) at eV/F = 0.205. The dashed and dot-dashed 
lines are for Pi{uj) and Pi{u), respectively. Panels (g) and (h) 
show plots of the DOS p(w) of the dashed line in panel (a) at 
eV = Ez (g), and at eV = —Ez (h). The other parameters are 
the same as for Fig.|S| 



the other hand, at negative bias voltages eV <—Ez, the 
spin-flip tunneling processes are between minority spin 
states. Hence the peak in the differential conductance is 
weak. 

A more precise explanation taking into account Kondo 
correlations can be given based on the spin-polarized 
nonequilibrium Kondo resonance. Figure Elg) shows 
p{uj) for eV = -Bz = 0.205 F. The two spin polarized 
Kondo resonances out of equilibrium split further into 
four peaks (solid line) by the energy eV Q. The lower 
two peaks at —Ez + eV/2 and —Ez — eV/2 are for spin- 
down (dot-dashed line) , related to the left and right lead 
chemical potential at ±6^/2, respectively. The peak at 
—Ez + eV/2 is larger than the peak at —Ez — eV/2, be- 
cause the spin-down electron in the QD can be screened 
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FIG. 8: The tunnel magnetoresistance TMR as function of 
the bias voltage V. Spin polarization factors are Pl = \ Pr \ = 
0.2. The further parameters are the same as for Fig. El 



by majority spins of the left lead with spin-up. The peaks 
at Ez ± eV/2 are for spin- up electrons (dashed line) , and 
a larger amplitude of the peak at Ez — eV/2 can be ex- 
plained in the same way as above. Under the condition 
eV — Ez, the positions of the two large peaks at Ez—eV/2 
and —Ez + eV/2 coincide with the two chemical poten- 
tials. Thus, tunneling electrons with energy between the 
left and right chemical potentials [shaded region in (g)] 
can profit from the larger DOS inside the QD. Conse- 
quently, the differential conductance Gap(V) at positive 
bias voltage eV = Ez is enhanced more strongly. Finally, 
Fig.Hh) shows the DOS p{lj) for eF= -i;z = -0.205 T. 
By repeating the same arguments as presented above, 
we can explain the fact that at negative bias voltage 
eV = —Ez the DOS p{lu) of the QD at energies between 
the left and right Fermi levels [shaded region in (h)] is 
lower, and consequently the peak of the differential con- 
ductance Gap{V) is smaller as well. 

The anomalous behavior of the nonlinear conductance 
determines the shape of the TMR. Figure |H1 presents the 
TMR defined as 



TMR = {Gp{V) - Gap{V))/Gap{V). 



(32) 



We find a large negative TMR at zero bias. This is in 
contrast to the TMR of ferromagnetic tunnel junctions, 
where one finds usually a maximum of the TMR at zero 
bias and then a decrease with increasing bias voltage 



D. Left-right asymmetry in system parameters 

In the previous subsection we showed that the spin- 
polarization of the Kondo resonance and the spin- 
accumulation for a nonequilibrium situation depend on 
the orientation of the lead magnetizations. We can ex- 
pect that a left-right asymmetry in coupling strengths, 
^Rt^^l, and spin polarization, Prj^Pl, or an asymme- 
try in the applied voltage could affect these properties 
as well and leave footprints on the nonlinear differential 
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FIG. 9: (a) The differential conductance for the P align- 
ment Gp{V) {Pl = Pr = 0.2). Solid lines, dashed lines and 
dotted lines are results for Fu/Fl = 1,0.5, and 0.25. (b) 
The differential conductance Gp{V) for rH/rL = 0.5 with the 
asymmetry factor ac — (dashed line) and 0.5 (dot-dashed 
line), (c) The differential conductance with the AP align- 
ment Gap{V) {Pl = —Pr — 0.2) and (d) the nonequilibrium 
local magnetization M. (e), (f) The total DOS p{u) in the 
AP alignment: (e) The solid (dashed) line is for r_R/r_L = 0.5 
at eV/r = -0.08 (0.07). (f) The solid (dashed) line is for 
rfl/ri=0.25 at eV/r = -0.13 (0.11). The other parameters 
are the same as for Fig. |S] 



conductance. In the following we discuss several exam- 
ples. 

Figure EI a) shows the nonlinear differential conduc- 
tance Gp{V) for the P alignment {Pl ^ Pr ^ 0.2) for 
varying values of the ratio of coupling strengths Tji/Tl, 
while keeping the total coupling F = Tl+Th unchanged. 
The values of the splittings are almost the same for all 
cases. The reason is the following: The spin splitting 
Eq. 1^ is expressed as the sum of the left and right con- 
tributions Ae = AeL + Ae/j (one can obtain the relation 
for Acr using Eq. ^ and adding a lead index r to the 
F and P). For P^ = Pjj ~ P, the value of the splitting 
Ae is proportional to P{Tl + F/j), which is independent 
of the ratio of coupling strength. The asymmetric ca- 
pacitive coupling does not alter the conductance G{V) 
by much. Figure Elb) shows curves for the asymmetry 
factor flc = (dashed line) and 0.5 (dot-dashed line) for 
^r/^l — 0-5 [the dashed line is the same as that in (a)]. 
The asymmetry factor just tilts the curve slightly. 

Figure |^c) presents the nonlinear conductance 
Gap{V) for the AP alignment for different values of 
Eii/Fi, and Oc = 0. For an asymmetric coupling, we 
find again a splitting of the zero bias anomaly. Such 
a splitting was observed experimentally |llj | (details are 
given in Sec.EJ- It can be understood because the spin 
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splitting due to the exchange interaction with the left 
Ae^ and right lead Ae^ cannot compensate each other 
AeL ^ — Aej^ for Tl ^ r^. Actually, with decreasing 
Ti^/Tl, Ae increases and thus the splitting is enhanced. 

We observe the asymmetry in the two peaks' amplitude 
similar to the case of the AP alignment with the Zeeman 
splitting [dashed line in Fig. [7|[a)]. Actually, this asym- 
metry can be explained almost in the same way. The only 
difference is that the large peak is at negative bias voltage 
[dashed and dotted Hues in Fig. EIc)] as opposed to the 
dashed line in Fig. [3a), because, Ae is negative whereas 
the Zeeman energy is positive. The solid (dashed) lines 
in panels (e) and (f ) show the total DOS p{uj) for the AP 
alignment at negative (positive) bias voltage. We can 
check that the larger DOS at negative bias voltage (solid 
lines) is responsible for higher conductance peak. 

Next, we will look at consequences of an asymmetry 
in left and right spin polarization factors Pl Pr- The 
dashed and dotted lines in Fig. ^| show the differential 
conductance G{V) (a) and local magnetization AI (b) 
for various spin polarization factors but for symmetric 
coupling strength F/j = F^. For each curve, the average 
value of polarization factors in the left and right leads 
was kept constant {Pl + Pr)/'^ = 0.1. The solid lines 
show the result for the equal polarizations {Pl = Pr — 
0.1). As the spin polarization of the left lead increases 
{Pl = 0.4 for dashed hues, Pl = 0.6 for dotted lines) 
the local magnetization increases (decreases) for positive 
(negative) bias voltage [Fig. Iior b)] and accordingly, the 
difference between the two peaks grows [Fig. llOf a')]. The 
solid (dashed) lines in panels (c) and (d) are the total 
DOS at negative (positive) bias voltage [(c) for Pl — 0.4 
and (d) for 0.6]. It confirms that with increasing Pl, 
the difference in the intensity increases, causing the large 
difference in peak heights. 



IV. DECOHERENCE EFFECT 

In this section, we will turn our attention to the ef- 
fects of decoherence. Since the transverse spin relaxation 
time, T2, is defined as the decay time of the off-diagonal 
components of the reduced density matrix P^ , one might 
identify it within the Markov approximation with the 
imaginary part of the corresponding off-diagonal compo- 
nent of the self-energy in Fourier space, 



I/T2 = -lina'^''^{e^~ea)- 



(33) 



The relevant frequency in the argument e^ — is the en- 
ergy difference between the states with spin a and the 
opposite spin a. From Eq. (|15|l . one can conclude, that 
I/T2 is closely related to the broadening of Kondo res- 
onance. When evaluating the self-energy we note that 
the lowest-order contribution to does not contain 
an imaginary part. In the following, we therefore, per- 
form the second order perturbation expansion in terms 
of F at zero temperature T = and era = £(7 — /ir < 0. 
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FIG. 10: (a) The nonlinear conductance G{V) and (b) the 
local magnetization M for various spin polarizations {Pl, Pr) 
and for symmetric coupling Fr — Tl- The average value of 
the polarization factors in the left and right leads are fixed as 
{Pl + Pr)/'2 = 0.1. The solid, dashed, and dotted lines are 
for Pl =0.1, 0.4, and 0.6, respectively. The total DOS p{uj) 
at el//r = -0.1 (solid line) and e\^/r = 0.1 (dashed line) for 
Pl =0.4 (c) and Pl =0.6 (d). The other parameters are the 
same as for Fig.|S| 



The first diagram of the second-order expansion, de- 
picted in Fig. |2Ib), is evaluated as 



(e^-Cff) = / duj'dhj" ] 

a' 

7+(c.")7.-(^0 

- e^' - cj' + cj" + zO+)(e<, - cj' + i0+) ' 



(34) 



For simplicity, we assume here a flat DOS with cutoff 
energy D, i.e., we replace the Lorentzian cutoff function 
by pc(w) e{ui + D)9{D - uj) with |e,.<,|<£). Then, the 
imaginary part of Eq. H34|) reads. 



— Im cr' 



{2^f 



■ 0{era - Erv) • (35) 



The second diagram of the second-order expansion is ob- 



tained from Eq. ^ as al%^'^~^\w) = ^-"""^ (-w)' 

In addition to the above two diagrams there are four 
other diagrams (Fig. Ill|) ; however their imaginary part 
vanishes for era < 0: 



,a.a(2-l)/ 



Im CT^'^' 



Im / duj'dw'^ 



j^{u;')j-{u;") 



t0^ 



CJ' - Eff + iO+ 



^0^ 



E 



F V I- 

^ TO r' a 

An 



{Cr'a) , Q{^rij) 



= 



(36) 

By adding the contributions from cr^'^^^ and crj'^*'^ 
one obtains the transverse spin relaxation rate in 2nd 
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FIG. 11: The four additional diagrams of second order expan- 
sion for the irreducible self-energy cr^]^. One bubble dresses 
each horizontal line. The other two diagrams, where one bub- 
ble is inside the other bubble are depicted in Fig.|3Ib). Here, 
dashed lines carry the energy — ^a- 



order. 



Lo-i Ra' l^La 



T, 



(2) 



(7(7' 



47r CLa ^R(7' 



(37) 



The rate is proportional to the sum of spin flip and spin 
preserving co-tunneling currents. This result is consistent 
with X^ffl/'^CT derived in Refs. ]^ ^'^'i-}^ defined in 
Eq. (5) of Ref. or Eq. (43) of Ref. p), where the 
decoherence effect due to the dissipative current flow for 
finite bias and/or magnetic field even at zero temperature 
was first pointed out. 

Because the RTA ignores crj^^, it predicts sharp con- 
ductance peaks split by the magnetic field 7] . As stressed 
in Refs. 5, 34], we cannot neglect the decoherence effect, 
since it could provide the cutoff energy for Kondo cor- 
relations and thus suppresses the Kondo effect. In sec- 
ond order perturbation theory, discussed above, the de- 
coherence effect appears to be weak, causing only slight 
changes of the DOS. However, a more thorough analy- 
sis reveals that second order perturbation theory is not 
enough: it does not account for Kondo correlations which 
also enhance the decoherence [3 |3^. In order to de- 
scribe this enhancement, we have to sum up infinite order 
diagrams again. 
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FIG. 12: (a) The off-diagonal component of the reduced prop- 
agator in spin space within NCA. Double lines on the Keldysh 
contour are determined by Dyson equations in (b). 

A convenient way to include diagrams ^up to the infi- 
nite order is provided by the NCA 2, 26,zf\. The naive 
application of NCA causes a spurious peak at the Fermi 
energy (see Sec. IJ), but the NCA for the off-diagonal 
component of the reduced propagator in spin space Tr?'^ 
is free from this ffaw as shown below. The diagram 
of 7T^[% within NCA is depicted in Fig. ^a). Double 
lines are 'full' propagators of the QD state and deter- 
mined by Dyson equations generating the infinite num- 
bers of 'non-crossing' diagrams [Fig. I12f b)]. We observe 
that diagrams in first and second order expansion of cr'^'g 
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FIG. 13: (a) The total and spin resolved DOS for paral- 
lel alignment {Pl = Pr = 0.2) within NCA for TrJ;^ (solid 
lines). The panel (b) is the magnification around the Kondo 
resonance. Dotted lines are the results without decoherence 
effect and correspond to lines in Fig.jSJc). Parameters are the 
same as those in Fig. 15] 



[Figs. Elb) and ^] are contained in the NCA diagrams. 
The expression is given by 



(J,cr / \ 



LU + iO+ — UJl + UJ2' 



(38) 



where Pxi^) is the spectral density of the full propaga- 
tor of the QD state \x), obtained by the following self- 
consistent equations: 



S<t(cj) 

Po{^) 
Eo(tj) 



1 

-Im- 



duj'%{uj - uj')po{uj'), 



1. 



-Im- 



1 



(39) 
(40) 
(41) 
(42) 



Figures IKSf a) and (b) show the total and spin re- 
solved DOS and their magnifications for parallel align- 
ment {Pl = Pr^0.2) within NCA for Tr^% (solid lines). 
By comparing the results without decoherence effect [dot- 
ted lines correspond to lines in Fig.E^c)], we can see to 
which extend the Kondo resonance is suppressed by the 
decoherence effect. 
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Though NCA for tt^'J describes the suppression of the 
Kondo resonance and also the spin spUtting, careful con- 
siderations of vertex corrections are necessary jsSj . More- 
over additional self-consistency conditions require exten- 
sive calculations. Even if analyzed numerically, this is 
beyond the scope of the present work, the emphasis of 
which is on spin splitting effects in magnetic systems. 
For further investigations for the decoherence problem, 
advanced approaches, such as the real-time renormaliza- 
tion group technique or the generalization of the 
perturbative renormalization group to nonequilibrium 
states '34^ may be appropriate. 



V. RELATION TO EXPERIMENTS 

In this section we will discuss the correspondence be- 
tween our calculation and very recent experimental re- 
sults. Pasupathy et al. studied the transport through 
a single Cgo molecule attached to ferromagnetic nickel 
electrodes in the Kondo regime. It was shown that the 
Kondo correlations exist even in the presence of ferro- 
magnetic leads. The zero-bias anomaly in the nonequi- 
librium conductance was split for P alignment of the lead 
magnetizations, in agreement with our calculations pre- 
sented in Figs.EIa), 7(a), 9(a), and 9(b). For AP align- 
ment the splitting of the zero-bias anomaly was absent 
in one case - similarly to Fig. [3c), and substantially re- 
duced in other cases, which can be interpreted as an effect 
of asymmetric coupling F^ ^ F^, similarly to Fig.|5Jc). 
The measurement of the nonequilibrium conductance for 
the parallel alignment for several temperatures demon- 
strate that the splitting of the Kondo resonance does not 
depend on temperature, in agreement with results pre- 
sented in Fig. EJa). Also the TMR signal showed very 
similar properties to those presented in Fig. |S| 

In a different set of experiments Nygard et al. re- 
cently observed in the transport through a carbon nan- 
otube, which acts as a QD, attached to normal leads, a 
zero bias anomaly at low temperature. This anomaly is 
split probably due to interaction effects between the car- 
bon nanotube and a magnetic catalyst particle. These 
results could be interpreted within the framework of our 
model as well. 



VI. SUMMARY 

In summary, we have analyzed the nonequilibrium 
Kondo effect in a QD coupled to ferromagnetic leads. 
We started form a systematic approach - the real-time 
diagrammatic technique. We proceeded in an extension 
of the 'resonant tunneling approximation' by introducing 
the self-energy of off-diagonal components of the reduced 
propagator in spin space, which was neglected in the orig- 
inal RTA. The self-energy was calculated to lowest order 
to describe spin-dependent quantum charge fluctuations. 
The approximation ensures spin and charge conservation. 



The local magnetization is obtained by solving the master 
equation accounting for quantum fluctuations. Our work 
gives a coherent description of spin splitting, spin accu- 
mulation, and Kondo correlations out of equilibrium. We 
calculated the nonlinear differential conductance for var- 
ious temperatures. Though the peaks hight of the split 
zero-bias anomaly is suppressed with increasing temper- 
atures, the spin splitting itself is robust against tempera- 
ture and leaves a clear structure in the nonlinear differen- 
tial conductance. Our approximation provides a reason- 
able description of the compensation of the spin splitting 
by an applied magnetic field. We also discussed asymme- 
tries in the zero-bias anomaly in the conductance split by 
the external magnetic field for the antiparallel alignment 
in term of the spin polarization of the Kondo resonances. 

We also discussed the effect of the asymmetry in 
system parameters. The asymmetry in the coupling 
strengths does not affect the spin splitting for the parallel 
alignment. For the antiparallel alignment, the spin split- 
ting is absent for symmetric set-ups. However, an asym- 
metry in coupling strength leads to a spin splitting be- 
cause the exchange interaction with the left and the right 
lead cannot compensate each other. The asymmetry is 
important especially for nonequilibrium experiments, be- 
cause it can induce a spin accumulation and affect the 
weights of the split Kondo resonance. We showed that 
the asymmetry in spin polarization factors of two leads 
can induce rather large asymmetries in the peak heights 
of the split zero-bias anomaly. 

Though the result looks qualitatively reasonable, our 
lowest-order perturbation theory for the self-energy falls 
short of a proper description of decoherence effects. We 
showed that the NCA for the off-diagonal component of 
the reduced propagator in spin space could account for 
the decoherence effect which suppresses the hight of the 
Kondo resonances. The further analysis of this prob- 
lem within the real-time approach remains a challenging 
problem for the future. We hope that the systematic 
formulation of nonequilibrium Kondo effects in magnetic 
systems presented above encourages further efforts to in- 
vestigate these effects in spintronics devices. 
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